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ABSTRACT 

We have used Aperture Masking Interferometry and Adaptive Optics (AO) at the Palomar 200" to 
obtain precise mass measurements of the binary M dwarf GJ 623. AO observations spread over 3 years 
combined with a decade of radial velocity measurements constrain all orbital parameters of the GJ 623 
binary system accurately enough to critically challenge the models. The dynamical masses measured 
are mi = 0.371 ± 0.015 M Q (4%) and m 2 = 0.115 ± 0.0023 M Q (2%) for the primary and the secondary 
respectively. Models are not consistent with color and mass, requiring very low metallicities. 
Subject headings: binary, stars: luminosity functions, mass functions, techniques: AO, interferometric 



1. INTRODUCTION 

The mass of a star, along with its metallicity and 
age, is the fundamental parameter that determines its 
position along an evolutionary track. Even if binarity, 
rotation, magnetic fields and other parameters also af- 
fect stellar interiors, what is known as the Vogt-Russel 
theorem remains an important rule, and the mass an 
essential parameter of stellar evolution. Binary sys- 
tems offer the ideal test to infer dynamical masses in- 
de pendent of the use of a stellar model such as those 
of iBaraffe etHI <H998), 

or an empirical mass-luminosity 
(hereafter M/L) relation. Such relations, available 
for both visible dHenrv et al.l Il999f ) and near infrared 
(jDelfosse et al.ll2000D are important astrophysical tools, 
fairly well constrained for intermediate mass stars. How- 
ever, the Solar neighborhood is domi nated by low mass 
stars, in both number and total mass (Henry' 1998), and 
as dust condenses in the atmospheres of these cool stars, 
the models meet new unknowns. There also are puzzles 
in the stellar structure of M dwarfs and the extrapolation 
of these models to substellar objects remains untested. 
Below 0.6 Mq, both M/L relations and models will ben- 
efit from model-independent determinations of high pre- 
cision stellar masses, which can be achieved by combin- 
ing radial velocime try measurements and high angular 
resolution imaging (Delfosse et al. 2004; Segransan et al. 
2000) . The complementarity of the two techniques yields 
substantial benefits, even in the regime where the time 
baselin e of the observations is s horter than the orbital 
period (jEisner fc Kulkarnill2002T ). 

We have observed known M dwarfs binaries with 
precise radia l velo city measurements published by 
iNidever et al.l (|2002h . These late type M dwarfs are ideal 
observing targets for imaging with the PALAO Adaptive 
Optics (AO) system and the PHARO infrared camera 
(|Hayward et aLlfeOOlh . optimized for the near infrared 
bands. At small angular separation (i.e. < 2X/d), the 
sensitivity of the detection of faint companions can be im- 
proved by combining AO with aperture maskin g interfer- 
ometry (fLlovd et al.fl2006l : iPravdo et al.ll2006h . This pa- 
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per provides the astrometry of the binary system GJ 623, 
successfully observed with both AO and AO + aperture 
masking. Combined with the radial velocity measure- 
ment, the astrometry provides precise dynamical masses 
(better than 2%) of the GJ 623 binary system. Combined 
with the J,H,K photometry, this measurement adds new 
constraints to the models and M/L relations. 

2. THE OBSERVATIONS 
2.1. The primary 

GJ 623 (aka LHS 417, HIP 80346) is a high proper 
motion M2.5 dwarf 8 pc from the Sun. It is a 
long known astrometric bina ry, first characterized by 
iLippincott fc B orgman ()1978l ) . G J 623 has proven to be 
a ideal test for diffe rent observational techniques, such as 
radia l velocimetry (M arcv fc Moorel [l989: Nidcve r" et alJ 
2002) and speckle interferometry ([McCarthy fc Henrvl 
Il987f ). This binary system has also been directly imaged 
in the visible with t he COSTAR corrected HST/FOC 
(jBarbieri et al.lll996h . 

Our observations of GJ 623 were performed with the 
PHARO instrument on the Palomar 200" telescope and 
with the NIRC2 instrument on the Keck II telescope. 
The companion of GJ 623 was detected in six observ- 
ing runs: September 2003, June 2004, January 2005 and 
February 2006, at Palomar, using J, H and K s filters, and 
June and August 2006 at Keck, with the 1.58 \xm H cont 
filter. Tables Q] and [2] respectively gather the astromet- 
ric and photometric measurements made at Palomar and 
Keck. 

Four of our observations consist of conventional direct 
imaging with AO. These images were dark substracted, 
flat fielded and analysed with a custom IDL program, 
using the latest PHARO and NIRC2 plate scale and 
orient ation characterization by iMetchev fc H illenbrand 
(|2004f h The location of the companion (angular separa- 
tion and position angle) and the constrast ratio are pre- 
cisely determined with a cross correlation of the images. 
Not surprisingly, doubling the diameter of the telescope 
(c./. the Keck data points) significantly improves the 
precision of the measurements and roughly reduces the 
error bars by a factor 3. 

The data also include observations in ape rture masking 
interf erometry with AO. Th is technique dTuthill et alJ 



I2000h . recently described bv iLlovd et ail (|2006h . is ap- 
propriate to the detection of faint companion at small 
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angular separation (typically less than 2 X/D), where 
direct AO imaging has so far proven difficult. 

Indeed, present AO systems focus on achieving high 
contrast at moderate angular separation > 4X/D. Below 
this limit, the variance of the speckle background dom- 
inates the photon n oise by several orders of magnitude 
(R acine et al][l9 99). This prevents us from calibrating 
the Point Spread Function (PSF), with a precision suf- 
ficient to discriminate faint companions from the bright 
star's residual speckes. 

This issue can be evaded by using intcrfcrometric tech- 
niques. By sampling a few spatial frequencies only, a 
mask located in the pupil plane permits us to com- 
pletely decompose the PSF into a finite set of Fourier 
components. Non-redundancy of the baselines passed 
by the mask ensures that each frequency is sampled 
only once a nd the visibilities can be used to form clo- 
sure phases ([Baldwin et alj|1986t iNakaiima et al.lll989l ; 
iReadhead et al.lll988f ). This observable rejects both at- 
mospheric noise and calibration errors of the wavefront 
sensor. The only drawback is the tranmission of such a 
mask: between 5 and 15 % for the ones used at Palomar. 
However, in the speckle noise limited regime, light loss 
does not result in a loss of sensitivity. 

Our analysis also uses a HST observation dating t o 
1994, and originally published bv iBarbieri et all (1996): 
an angular distance of 330 ± 20 mas and a position angle 
of 7.0±2.6° on June 11, 1994. Even though this is inferior 
to the precision achieved with Palomar and Keck, which 
are much larger telescopes, the fact that this observation 
was performed more than 10 years ago adds an important 
constraint on the period of the binary. 



Julian Date 


Band 


Pupil 


Sep. 


PA 


(-2450000) 




(see text) 


(mas) 


(deg) 


2896.6 


Ks 


9 hole 


240.4 ± 9.7 


79.31 ± 2.0 


3163.8 


[Fe II] 


18 hole 


340.5 ± 20.4 


49.16 ± 4.3 


3402.1 


K s 


full ap. 


350.7 ± 2.2 


28.66 ± 0.5 


3780.0 


J,H,K S 


full ap. 


257.1 ± 3.1 


-12.76 ± 0.6 


3909.5 


Hcont 


Nirc2 


176.7 ± 1.1 


318.2 ± 0.2 


3962.3 


Hcont 


Nirc2 


138.7 ± 0.9 


295.4 ± 0.2 



Table 1 

astrometric measurements at palomar and keck: 
angular separation and position angle of gj 623 b. 



The Hipparcos parallax it = 124 . 34 ± 1.16 mas origi- 
nally published bv iPerrvman et ail (119971) was based on 
the ra dial velocimetry measurement of iMarcv fe Moore 
(1989). The parallax was recalculated bv I Jancart et al. 
(2005) after the publication of an improved radial veloc- 
ity curve bv lNidever et alJ (|2002l ). It is this revised par- 
allax 7r = 125.81 ± 1.19 mas (Pourbaix, D., priv. comm.) 
that we adopt here. The following absolute magnitudes 
can be deduced for GJ 623 AB: Mj = 7.137 ± 0.052, 
M H = 6.640 ± 0.051 and M K = 6.414 ± 0.052. 

3. METHOD: EXTRACTING ORBITAL PARAMETERS 

The starting point in the determination of the orbital 
parameters of the GJ 623 system, is the set of radial ve - 
locity measurements published by iNidever et alJ (|2002f ) 
and available in the online version of the Ninth Cat- 
alogu e of Spectroscopic Binary Orbits (jPourbaix et al.l 
l2004f ) at |http : //sb9 . astro . ulb . ac . be j The parame- 
ters deduced from the Keplerian fit to the data are sum- 
marized in table [51 They have been directly used by 
I Jancart et all (|2005l ) to re process the Hipparcos Interme - 
diate Astrometric Data of Ivan Leeuwen fe Evans! (|1998l ) 
for a complete characterization of the dynamical ele- 
ments of the system. 

The approach detailed in this paper consists of com- 
bining together, with no a priori assumptions, the orginal 
radial velocity data and our astrometric observations. 

3.1. Keplerian Orbits 

One uses the standard 2-body solution to parametrize 
the location of the companion orbiting GJ 623 A. To an 
observing date t, one associates an angle M called the 
mean anomaly: 

M(t) = T (t-T P ), (1) 

where P and Tp respectively represent the orbital period 
and the epoch at the periastron passage. In the orbital 
plane, the (x, y) coordinates can simply be expressed as 
a function of another angle, the eccentric anomaly E. It 
is the angle between the direction of the periastron and 
the current position of the companion, projected onto 
the ellipse's circumscribing circle perpendicularly to the 
major-axis, measured at the center of the ellipse (c.f. fig. 



Julian Date 


Filter 


Amag 


(-2450000) 






2896.6 


K s 


3.051 ± 0.826 


3163.8 


Fe II 


2.903 ± 0.476 


3402.1 


K s 


2.604 ± 0.047 


3780.0 


J 


2.691 ± 0.038 




H 


2.860 ± 0.039 




K s 


2.789 ± 0.014 


3909.5 




2.794 ± 0.033 


3962.3 




2.781 ± 0.016 



Table 2 

Photometric measurements of GJ 623 B. 



From the 2MASS catalog, the apparent magnitudes of 
GJ 623 AB are J = 6.638 ±0.024, H = 6.141 ±0.021 and 
K = 5.915 ±0.023. 



M(E)=E-esmE (2) 
x{E)=a(cosE - e) (3) 

y{E) = ayfl - e 2 sin£, (4) 

where a and e respectively represent the semi-major axis 
and the eccentricity of the orbit. Eq. O is called Kepler's 
equation. One solves this equation i.e. finds the value 
of E associated with a given M by using the following 
classical iterative algorithm: 

E = M 

E n+l =M + esmE n . (5) 

3.1.1. Radial Velocity Orbital Models 

One can use the standard 2-bodies solution (eq. [3] and 
|4|) to compute the coordinates of the primary component 
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is the angle between the direction of the periastron and 
the current position of an object on its orbit (P on fig. 
[J), measured at the focus of the ellipse. 

Fig. [T] illustrates the one-to-one correspondance be- 
tween eccentric and true anomaly. The projection of the 
point P on the major axis provides the following relation: 



r cos v = a(cos E — e) . 



(11) 



Together, eq. [10] and [IT] lead to the primary to sec- 
ondary distance r and the true anomaly v, as functions 
of the eccentric anomaly E: 




Fig. 1. — Parametrization of an elliptic orbit. Geometric relation 
between the excentric anomaly E and the true anomaly v. 



in the orbital plane, relative to the center of mass of the 
system: 



x P (E) 

y P (E) 



m 2 
rriT 



x(E) 

y(E) 



(6) 



where mj- and m 2 respectively stand for the total mass 
(primary + secondary) and the mass of the secondary. In 
the orbital plane, the velocity vector can be calculated 
with the following partial derivative: 
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dE 8M 

x m x lif (7) 



This derivation leads to the following x and y compo- 
nents of the velocity: 



v x {E) 

Vy(E) 



m 2 



2na 



m,T P(l — ecosE) 



sinE 
- e 2 cos E 



(8) 

The Radial Velocity (hereafter RV) is the component of 
this velocity projected on the line of sight. Its expression 
therefore involves both the argument of the periastron 
luq and the inclination of the system i. The RV is not 
sensitive to the orientation of the system on the sky, i. e. 
to the value of the argument of the ascending node Qi: 

RV(E,ujq, i) = [v x (E) sinojo + v y (E) cosoj ] sini + V0) 
where Vq is a constant term, the mean RV. 

3.1.2. Astrometric Orbital Models 

Unlike the case of RV data for which equations are 
expressed from the center of mass of the system, the 
astrometric measurements gathered in Table Q] provide 
the instant position of the companion relative to the pri- 
mary. In this frame, the trajectory of the companion is 
an ellipse, whose focus is the primary component. The 
parametric equation of this ellipse is: 



r{v) 



t (W) 



1 



ecosv 



(10) 



where a is the semi-major axis, e, the eccentricity (0 < 
e < 1) and v, an angle called the true anomaly, which 



r = a(l — ecosi?), 
cos E — e 



1 — e cos E 



(12) 
(13) 



At a given observing date t, one needs once more to 
solve Kepler's equation (c.f. eq. ^ to determine the 
corresponding eccentric anomaly E. The location of the 
companion along the orbit is provided by eq. [12] and [T3l 

Contrary to the radial velocity, which is the component 
of the velocity projected on the line of sight, one mea- 
sures here the position of the companion projected on 
the celestial sphere. Right ascension a and declination S 
of the secondary (relative to the primary) are given by 
the following relations: 



a = r x [cos(v + ujq) sin fij + 
sin(^ + ujq) cos i cos fii] 

6 = r x [cos(z/ + uiq) cos Qi — 
sin(^ + ujq) cos i sin . 



(14) 
(15) 



3.2. x 2 fitting 

We use the model presented in the previous section, to 
fit a 9 parameter model to 26 observables. The observ- 
ables are 2 coordinates for 7 astrometric data points and 
12 radial velocities. The 9 parameters are the 6 orbital 
elements, i.e. the semi-major axis a, the eccentricity e, 
the longitude of the ascending node u>o, the inclination i, 
the argument of periastron f^i and the orbital period P, 
plus the RV offset Vq and the semi-amplitude of the RV 
curve K\. 

In the case of a conventional analysis of RV data, the 
main observable is the semi-amplitude K±, which once 
combined with the period and the eccentricity of the or- 
bit, provides the mass function: 



f(M u M 2 ,i) 



M| sin i 3 
(M T ) 2 



(16) 



This combined analysis RV+astrometry allows us to 
separate the geometrical effects from the semi-amplitude. 
This produces another composite observable, somewhat 
simpler than K\ , a "pseudo-amplitude" , whose formal 
expression is: 



M 2 2i:a 



(17) 



This is the parameter that will be used to determine the 
dynamical masses of the two components of the binary. 
The final result of this 9 parameter fit, with 17 degrees of 
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GJ 623 Radial Velocity curve 
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Fig. 2. — Orbit of GJ 623B. Measurements and associated uncer- 
tainties are represented in blue. Points 1- 4: PHARO o bservations, 
point 5: Hubble observations reported in {Barbicri ct al. 199j|) and 
points 6-7: Nirc2 observations. In green, the Keplcrian fit to the 
orbit. The red line is the line of nodes, the point P marks the 
position of the Periastron. 



freedom is represented on Figure[5]for the astrometry and 
[3] for the radial velocimetry. Our solution exhibits a final 
reduced \t — 1-03, which is, despite the heterogeneity of 
the data, close enough to unity to ensure confidence in 
our estimation of the error bars. It is dominated by the 
velocimetry, despite having fewer measurements than the 
astrometry. 

The confidence interval of each of the parameters is de- 
termined by analysis of the likelihood function. If we as- 
sume that the the noise associated to our measurements 
is Gaussian, and that our parameters are independent, 
the likelihood can be approximated by: 

( x 2 \ 

L (parameters) oc expl —J. (18) 

The computation of a 9-dimcnsional likelihood func- 
tion requires a lot of CPU time. This difficulty may be 
circumvented by confining the search to a subset of the 
total space. This is achieved by fixing the values of cer- 
tain parameters and calculating the joined likelihood for 
the remaining parameters. However, one needs to check 
a posteriori that the validity of the hypothesis of inde- 
pendent parameters. 

This analysis reveals that among the 9 parameters, 
only the argument of the periastron luq and the longi- 
tude of the ascending node fii exhibit significant cor- 
relation. These two parameters are constrained within 
0.5° (c.f. Table [3]). The hypothesis of independent pa- 
rameters is valid, as expected from our good coverage 
in both radial velocimetry and astrometry. The uncer- 
tainty associated to each parameter is taken equal to the 
standard deviation of its associated likelihood function. 



Fig. 3. — Radial Velocity curve of GJ 623A. The mea surements 
and associted uncertainties reported by iNidever et alj|2002l) are 
represented in blue. The green solid line represents the Keplerian 
orbit with the parameters derived by \ 2 fitting. 



The Gaussian-like likelihoods of the parameters used in 
the determination of the dynamical masses are shown on 

Fig. rj 

4. CHARACTERISTICS OF GJ 623 AB 

The orbital parameters and the confidence intervals 
derived from the likelihood analysis are summarized in 
Table [3] This combined analysis significantly improved 
the constraints on t he epoch at periastr on passage, from 
T P = 2451298±10 (jNidever et al.ll2002j ) to the new value 
T P = 2451313.3 ± 0.6. The RMS of the Radial Velocity 
fit is 56 m/s for a semi-amplitude of K\ = 2.01 ± 0.01 
km/s, which is comparable to the original fit (RMS=51 
m/s for K x = 2.08 ± 0.04 km/s). 



parameter 


Nidever 


This work 


n (mas) 




125.81±1.19 


a (AU) 
a (mas) 




1.894±0.019 
237.28±0.88 


e 


0.67±0.01 


0.631±0.002 


i (dcg) 




154.0±0.1° 


fl! 




98.5±0.47° 


UJo 


251±1° 


248.68±0.46° 


P (days) 


1366.1±0.4 


1365.6±0.3 


T P (reduced JD) 


1298±10 


1313.3±0.6 


V (km/s) 


-27.654±0.3 


-27.729±0.005 


A 1 (km/s) 




3.57±0.01 



Table 3 
Orbital Elements 



4.1. Dynamical masses 

From the data summarized in Table [31 one can deter- 
mine the dynamical masses of both components of the 
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Fig. 4. — Likelihood functions of semi-major axis a, period P and radial velocity pseudo-amplitude A\, used to determine the dynamical 
masses of the GJ 623 system. 



system. The revised parallax figuring in the table takes 
the binarity into account. The semi-major axis a, ex- 
pressed in AU, once combined with the orbital period, 
gives the total mass Mt and the associated uncertainty 

0~T- 



M T = a 3 /P 2 



(19) 

ctt/M t = v/(3a a /a) 2 + {2a P /Pf. (20) 

The likelihood analysis also determines the "pseudo- 
amplitude" defined by equation [17] and its associated un- 
certainty. 

This composite parameter, combined with the total 
mass, allows independent determination of the mass of 
the secondary M^: 



Mi = M T x 



Ax x P _ Ax x a 2 
2ira ~ 2nP 



(21) 



2 



/M 2 = ^(a Al /A 1 r + {op IP) 2 + {2a a /a)\ (22) 



and the mass ratio Mz/Mt" 



M R =^ 
M T 



Ax* P 
2ixa 



/M R =^(a Al /Ax) 2 + (ap/P) 2 + {ajaf 



(23) 
(24) 



A similar analysis is possible for M\. Errors on the 
masses are dominated by the uncertainty on the semi- 
major axis, itself dominated by the error on the Hippar- 
cos parallax. As a consequence, with a fractional error of 
1 %, the mass ratio is better constrained than the mass 
of the secondary. The dynamical masses of the GJ 623 
system are summarized in Table |U 



Quantity 


Value 


Total Mass 


M T = 0.486 ± 0.015 M Q 


Primary Mass 


Mi = 0.371 ± 0.015 M Q 


Secondary Mass 


M 2 = 0.115 ± 0.0023 M Q 



Table 4 
Dynamical Masses 



4.2. Color, Metallicity & Kinematics 

The multi- wavelength contrast ratios found in Table [2] 
may be used to decomposed the observed combined mag- 
nitude of the binary system into magnitudes for individ- 
ual components. One determines the following absolute 
magnitudes: Mj = 7.224 ± 0.052, M H = 6.719 ± 0.051 
and M K = 6.495 ±0.052 for GJ 623 A and Mj = 9.915 ± 
0.065, M H = 9.512 ± 0.052 and M K = 9.269 ± 0.053 for 
GJ 623 B. 

With J~K color indices of 0.729 ± 0.074 and 0.646 ± 
0.084 respectively for the primary and the secondary, the 
GJ 623 AB system is bluer than the lDelfosse et al.1 (|2000l ) 
empirical M/L relation. Therefore, we suspect GJ 623 
is of low metallicity. From the Hipparcos proper mo- 
tion (1145.38,-452.37) mas/yr (jPerrvman et alj I1997D 
and the RV offset Vb = —27.4 km/s determined from 
our model fit, one can calculate the Galactic space ve- 
locity (U, V, W) = (-33, 14, -41) km/s after correction 
for standard solar motion. This velocity is consis- 
tent with an ol d disk population, an d therefore of subso- 
lar metallicity (Chiba & Beers 2000), which is consistent 
with both components being bluer than the average field 
object. 

Figure O compares the location of both components 
of the GJ 623 system in a mass-lumi nosity diagram to 
the lo w mass population II models of iMontalban et all 
(2000) for different metallicities. The trend we observe 
with these models (c.f. Fig. [5]) supports GJ 623 AB bee- 
ing of subsolar metallicity. The model that best matches 
our measurements for the primary is for [M/H] = —1.0. 
For the secondary, the best model predicts a slightly 
lower mass of 0.110 ± 0.001 M Q for [M/H] = -2.0. This 
large discrepancy in metallicity is inconsistent with the 
assumption of a co-eval binary. The very low metallic- 
ity for the secondary would be consistent with GJ 623 B 
belonging to the galactic halo, which is unlikely accord- 
ing to the kinematics. Therefore, we conclude that the 
models do not adequately fit the data. 

5. CONCLUSION 



4 The sign convention is the one of the IDL astrolib gaLuvw pro- 
cedure. Note that the literature is confusing on this matter. For 
instance, McCarthy & Henrv ( 1987) provide numbers of compara- 
ble magnitude but with opposite sign for all (U,V,W) components. 



6 



F. Martinache et al. 



« 8 



10 



11 





[M/H] = -1.0 : 
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Fig. 5. — Mass-luminosi t y dia gram for the GJ 623 system in J, H and K bands. Comparison with the low mass population II stars 
models by Montalban ct al. (20Q3) in the 0.1-0.4 solar mass range, for different metallicities. 



The observation of binary systems is the only way to 
measure unbiased masses. As shown in this paper, com- 
bined with excellent radial velocity measurements, a few 
high angular resolution images provide sufficient infor- 
mation to constrain the range of possible masses below 
the 5% precision that is required to seriously challenge 
the models at the low end of the main sequence. 

The application of precision RV methods developed for 
planet searches, combined with AO will provide a com- 
plete characterization of the stellar structure of the lower 
main sequence. 
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